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Abstract 

We present Hall effect measurements on highly oriented pyrolytic graphite that indi- 
cate the occurrence of the integer quantum-Hall-effect. The evidence is given by the 
observation of regular plateau-like structures in the field dependence of the trans- 
verse conductivity obtained in van der Pauw configuration. Measurements with the 
Corbino-disk configuration support this result and indicate that the quasi-linear and 
non-saturating longitudinal magnetoresistance in graphite is governed by the Hall 
effect in agreement with a recent theoretical model for disordered semiconductors. 
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Recent experimental[l,2,3,4,5,6] and theoretical [7,8,9] work demonstrates the 
occurrence of correlated phenomena in the electron system of graphite, indi- 
cating that the classical view of the physics of its magnetotransport properties 
is not adequate. Besides, the quantum Hall effect (QHE) [1] as well as evidence 
for Dirac fermions (massless particles due to the linear dispersion relation) [5] 
have been reported for bulk samples of highly oriented pyrolytic graphite of 
high quality. Recently published results obtained on bulk graphite [10,11], a 
few layers thick graphite [12], and in graphene (single graphite layer) [13,14] 
have confirmed both observations and showed that in graphene the QHE is 
anomalous in that the contribution of the Dirac fermions makes the plateaus 
to occur at half-integer filling factors. Certainly, the understanding of the 
electronic properties of graphite is necessary for the development of advanced 
technology of the graphite-related nanomaterials. 
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The experimental work we report here is based partially on data presented in 
an international meeting [15]. These data clearly demonstrate the quantization 
of the Hall effect and in particular the existence of an even QHE in graphite, 
i.e. a transverse voltage that is symmetric on magnetic field reversal. The 
occurrence of the QHE in bulk graphite indicates that the couphng between 
graphene layers is much smaller than the one assumed by the Slonczewski- 
Weiss-McClure theory [16] but perfectly agrees with the results obtained by 
Haering and Wallace [17]. Also, the experimental results we present here in- 
dicate a possible solution to the longstanding problem of the quasi-hnear and 
non-saturating magnetoresistance (MR) [18] providing evidence that the MR 
at high enough fields is governed by the Hall resistance, in agreement with a 
recent theoretical model [19]. 

The occurrence of the QHE is not restricted to perfectly two-dimensional 
(2D) systems, as measurements in highly anisotropic layered systems [20], or- 
ganic Bechgaard salts [21], Molybdenum bronzes [22], layered crystals of the 
type (Bio.25Sbo.75)2T3 [23], (TMTSF)2AsF6 system [24], as well as in 200-layer 
quantum-well structures [25] indicate. Nevertheless, it is surprising that for 
the paradigm of strongly anisotropic systems, graphite, no clear evidence for 
the QHE has been published till 2003 [1,2,3]. The main obstacle to obtain 
clear evidence for the QHE in graphite is the sample quality, i.e. the internal 
short circuits between graphene planes caused by defects and impurities. As 
a characterization of the sample quality and orientation of the samples one 
may use the full width at half maximum (FWHM) of the rocking curves. In 
general we observed that the smaller the FWHM and the distance between 
voltage electrodes, the clearer are the plateaus in the Hall effect. The sam- 
ples measured in this work are highly oriented pyrolytic graphite (HOPG) 
obtained from Advanced Ceramics (AC sample) and from the research insti- 
tute "Graphite" in Moscow (sample HOPG-3), with FWHM equal to 0.40° 
and 0.64°, respectively. 

Angle dependent magnetoresistance (MR) oscillations as well as a ratio of 
the order of 10"^ between c-axis and in-plane resistivities in highly oriented 
graphite samples [2] indicate a much larger anisotropy, a much smaller coupling 
between the graphene layers and a greater infiuence of the two- dimensionality 
of the sample on the transport properties than earlier studies have reported. 
Among these findings is the magnetic-field induced metal-insulator-transition 
(MIT), which shows significant similarities to the one discovered in dilute 2D 
electron systems [26]. As in graphite [1], in some of these 2D systems the 
MITs have been shown to be connected to quantum-Hall-insulator-transitions 
at high magnetic fields [27]. After these experimental observations and the 
first theoretical prediction for the QHE in graphite [28], which was followed 
by further developed theoretical studies [29,30], the occurrence of quantum- 
Hall-states in graphite was already expected. 
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The longitudinal resistance as a function of magnetic field of all samples has 
been measured with the conventional ac-method, where four electrodes made 
of silver paint are placed on top of a surface parallel to the graphene layers. To 
calculate from this measurement the absolute resistivity (longitudinal as well 
as transverse) it is necessary to know the effective penetration depth of the 
current A. Independent measurements at room temperature using electrodes 
at different positions of the same sample reveal that for macroscopic, bulk 
samples A is in the range 1 /xm< A < 10 yum. With a conventional electrode 
geometry a good estimate of the resistivity is given by p ~ bW/Id where 
b is the width of the sample, / the applied electrical current, and V and d 
are the voltage and distance between the inner contacts. However, part of 
the data presented here were obtained using a van der Pauw configuration 
consisting of four point- like electrodes (1 ... 4) placed at the corners of the 
main sample surface. In the van der Pauw configuration the typical surface 
between contacts fixed at the corners of a parallelepiped was ~ 10 mm^. 
If we apply a current lu between contacts 1 and 2 and we measured the 
voltage V34 between contacts 3 and 4, the longitudinal resistivity is given by 
the equation p^x = 4.53AV34//12 [31]. The transverse resistivity is obtained 
from this arrangement measuring the voltages V42 and Vis with the input 
currents /13 and I42, respectively, using p^y — 0.5A(V^2/-^i3 — 143/-^42) [31]- The 
longitudinal and transverse conductivities axx and axy are calculated from the 
measured resistivities inverting the two dimensional matrix, i.e. 



Pxy 



Pxx Pxy 



'xy 



Pxx Pxy 



The magnetic field was always applied perpendicular to the graphene layers 
(the main sample surfaces). For all geometries the applied current was of the 
order of 1 mA, which is in the ohmic regime and for which no heating effects 
are observed. The resistance was measured with a LR-700 resistance bridge 
from Linear Research Inc., which works at low ac- frequency. 

The essential signature of the QHE is the occurrence of plateaus in the Hall 
conductivity axy{H) or resistivity pxy{H) isotherms. The transverse conduc- 
tivity value at the position of the apparent first plateau-like feature at H 
5T/po is for this sample (HOPG-3) [1,2,3] axy{H ~ 5T/po,T < 2 K) ~ 
0.27/A r2~^m^^, where A is given in meters. Assuming an exponential decay 
of the current with the distance from the surface and taking into account that 
the current penetration depth A >> c/2, the distance between the graphene 
layers c/2 — 3.35 A, the contribution of a single graphene layer can be written 
as 
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If we assume that the first plateau occTirs when = axy{c/2) = e^/h — 
3.86 X 10~^ using our data we estimate a penetration depth A = 2.35 fim, 
which is within the expected range. Figure 1 shows both conductivities for the 
HOPG-3 sample obtained from the measured resistivities and calculated using 
A = 2.35 iim. This figure shows that regular plateaus are clearly observable 
in the a^y{H) (or in axy{H)) measured in the van dcr Pauw configuration 
and plotted as a function of the inverse applied field. From the data in Fig.l 
we obtain that the transverse conductivity for a single layer at the apparent 
first plateau agrees within experimental error (±50% due to the uncertainty 
in the absolute value of A) with the separation between the plateaus Acr^ ~ 

3 X 10-5 n-^. 

Considering the simplicity of the model and the geometrical errors involved in 
the determination of the absolute values the quantitative agreement with the 
expected value for an integer-hke QHE seems reasonable. Within experimental 
error, the value obtained for (i^y depends on sample and also on the electrode 
arrangement, indicating that the disorder and current distribution may have 
an infiuence. From the sign of the Hall signal we conclude that the carriers 
have the same electronic charge as the electrons. 

A vanishing longitudinal resistivity Pxx{H) in the vicinity of the fields where 
the Landau levels are filled is a concomitant of the QHE, which appears in 
2D systems at very low temperatures [32]. However, p^xiH) (or a,j.x{H)) of 
graphite does not show any minima or a decreasing behavior with field, see 
Fig.l. The positive MR of graphite is very large in clear contrast to the rela- 
tively small MR in most of the 2D systems showing the QHE. Moreover, it is 
nearly linear in field without any sign of saturation (see Fig.l). High resolu- 
tion measurements of the MR of HOPG samples indicate a quasi-linear field 
dependence at fields as low as ~ 10"^ T with an anomaly in the field expo- 
nent at the MIT at HqH ~ 0.1 T [33]. Note that the magnetoconductance at 
high fields, where the Shubnikov - de Haas (SdH) oscillations due to the Lan- 
dau level quantization occur, shows a striking similarity with the transverse 
conductivity, see Fig. 1. 

The origin for the linear MR (LMR) in graphite as well as in other semimetals 
is a matter of current discussion [18,19,34]. The simple two band model [35] 
does not explain the quasi-LMR observed at low and high fields, unless one as- 
sumes ad-hoc field dependences for the electron and hole mobility. Abrikosov 
proposed that the LMR takes place in the Landau level quantization regime 
and above the field Hql that pulls carriers into the lowest Landau level. In 
order to account for the low values of the characteristic field {Hql ~ 10"^ T), 
the graphite-like energy spectrum with a linear dependence of energy on mo- 
mentum has been assumed [18,34]. The theory clearly states that LMR should 
be a generic property of graphene, which is the model system. The evidence we 
provide for an integer-hke QHE casts doubts whether the transport properties 
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Fig. 1. The components of the conductivity tensor per graphene layer vs. inverse 
field for the sample HOPG-3 measured in van der Pauw-configuration at 0.1 K (o) 
and 2 K (•). The conductivities were calculated using Eqs. (1) and (2) and assuming 
a penetration depth A = 2.35 nm. The absolute values of the conductivities axx and 
(T xy can be easily obtained dividing the plotted values by 3.35 x 10 ^'^m. Smaller 
(larger) values of the penetration depth A increase (decrease) the absolute values of 
the conductivities. In the off-diagonal clement regular plateaus are seen, indicative 
of the quantum-Hall-effect. Wc note that in this configuration the transverse Hall 
voltage changes sign when the field is reversed. 

of graphite can be appropriately described by this "quantum magnetoresis- 
tance" model added to the fact that a quasi-linear MR is observed also at low 
fields. 

Experimental evidence based on a large number of transport measurements in 
different HOPG samples indicates that the internal disorder of a graphite sam- 
ple influences the field and temperature dependence of the transport properties 
in such an extent that the QHE can in some cases transform in a linear field 
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dependent Hall effect with minima at specific fields [1]. Moreover, different 
contact distributions on the same sample show also a small but striking influ- 
ence on the transport properties. Therefore these samples should be considered 
as disordered semimetals. Regarding the quasi-linear and non-saturating MR 
in disordered semiconductors, Parish and Littlewood showed recently that this 
behavior might be provided internally by the Hall effect [19]. In what follows 
we demonstrate experimentally that the Hall-like signal can be measured in 
a graphite sample at large enough fields in a voltage contact configuration 
where for the homogeneous case this signal is not expected. 

The sample for a Corbino-disk experiment has been prepared by means of 
e-beam lithography using a 200 nm gold layer sputtered on the sample (AC) 
main surface, as shown in the inset of Fig. 2. In the Corbino disk geometry 
the diameter of the outer ring was 1.4 mm and 0.1 mm its width as well as 
the diameter of the electrodes. In this configuration, the current is applied 
between the inner point-like electrode and four contacts along the outer ring- 
like electrode to assure homogenous current input. The voltage is measured 
between any of the other contacts. In a homogenous, isotropic (along the 
planes) 2D sample, the application of a constant current / would result in 
a two-dimensional radial current density jr — //27rr at a point located at 
distance r from the center, which is deflected by the magnetic fleld by the Hall- 
angle 6 with tan(^^) = Cxy/cxx- The result is a current path along logarithmic 
spirals, consisting of the radial component jr and an azimuthal (circulating) 
component given by [36] = \jrO'xy/o'xx\, see the inset in Fig. 2. If we neglect 
any influence of the Hall component, the voltage between any two contacts at 
the distances ri and r2 from the center is then given by 



Here an is the longitudinal in-plane conductivity. As expected, the voltage 
vanishes between contacts located at the same distance from the center. The 
circulating current component is caused by the Lorentz force perpendicular to 
the electric fleld and therefore it is a non-dissipative current, which does not 
contribute to the voltage. The aim of the experiment is to show that when 
for a given pair of electrodes inside the Corbino disk the "longitudinal" -like 
signal is minimized (i.e. ri ~ r2, see Eq. (3)), still a voltage is measured and 
this resembles the Hall signal. 

Figure 2 shows the fleld dependence of the inverse conductivity obtained from 
voltage measurements between various contact pairs. Two of them (A-C and 
G-F) provide the longitudinal magnetoconductivity and three pairs at nom- 
inally the same distance from center (B-C, E-F and C-D) should show no 
voltage. In contrast to the expectations for a homogeneous sample, the volt- 
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Fig. 2. Inverse diagonal component of the conductivity tensor vs. applied field for 
the AC sample measured in Corbino-configuration at 2 K. The conductance g ~ a\j 
is calculated from Eq.(3) assuming a relation of the distances of the voltage con- 
tacts X = r^lrx which leads to the best scaling onto the curve for the contact pair 
A-C (see below). The current is applied between the inner electrode and the outer 
ring; voltage is measured between pairs of the other electrodes. The curves were 
obtained from the voltage contacts: (+) A-C, x = 1.5; (X) G-F, x = 1.957; (□) 
B-C, X = 1.111; (O) E-F, x = 0.915; (•) C-D, x = 1.009. The electrode designations 
correspond to the ones in the inset. Inset: View from top onto the Corbino- configu- 
ration with designations of the electrodes as used in the text. The current path along 
a logarithmic spiral is schematically shown together with the radial component jr, 
the azimuthal component and the Hall-angle 9. 

age between pairs located at nominally the same distance from the center is 
not negligible. The first explanation for this observation would be that this is 
due to a misplacement of the contacts giving rise to a radial distance between 
them. To check this, the data in Fig. 2 is presented as follows. For the contact 
pair A-C (see inset in Fig. 2) the conductance (7 ~ an is calculated from (3) 
assuming the nominal value for x = r2/ri = 0.45 mm/0.30 mm=1.5. For the 
other pairs, x is chosen in such a way, that the curves scale onto the one for 
the pair A-C. 

As can be seen in Fig. 2 the curve for the pair G-F scales very well assuming 
X = 1.957, corresponding to a shift of the inner contact from the nominal 
one by 70 /xm, caused by the limited accuracy of the preparation (nominal 
electrode diameter is 100 /xm). The same holds for the pairs B-C and E-F, 
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Fig. 3. Data obtained from the contact pair C-D in the Corbino-configuration. Top: 
Curve taken at 2 K. Bottom: Curve taken at 5 K at positive (•) and negative (o) field 
directions. Note that the voltage does not change sign when the field is reversed. 

and therefore we conclude that those signals are mainly determined by the 
"nominal" longitudinal MR. The pair C-D, however, has the smallest nominal 
misplacement {x = 1.009), but shows the most significant deviation from the 
assumable purely longitudinal MR curve A-C. Interestingly, C-D is the contact 
pair where we observe the clearest platcau-likc structures similar to the ones 
measured in the van der Pauw configuration, as shown in the upper part of 
Fig. 3. The periodicity of the plateaus agrees very well with that obtained 
from the van der Pauw configuration. The inverse-field period of the plateaus 
is A{1/ fioH) ~ 0.17 T~^. Using the Onsager relation [37] for a 2D electron 
system A{1/ fioH) = e/'Khn2D we obtain a 2D electron density for our HOPG 
samples n2D — 2.85 x 10^^ m~^. This value is comparable with literature values 
for the 3D electron density 713^1 ~ 3 x 10^^ m~^ ~ n2D2/c. Prom this evidence 
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we conclude therefore, that the measured voltage between the electrodes C-D 
at large enough fields is not related to the longitudinal MR. 

The origin for the unexpected Hall signal at the pair C-D is compatible with 
the disordered semiconductor model of Ref.[19]. In this model the LMR is 
governed by the Hall resistance of the current paths. Note that due to the 
influence of pxy in Pxx)<^xy ~ ^/ Pxy However, this signal is independent of 
the field direction (switching the field direction will also switch the transverse 
direction of the current within the disordered network). That means that the 
measured MR at C-D should be field direction independent in case the Hall 
current paths are perfectly symmetric against field inversion. To check this 
we performed measurements at the two field directions, i.e. +z and —z. The 
results shown in Fig. 2 and in Fig. 3 indicate that the observed voltage does 
not change sign and does depend only slightly on the field direction probably 
due to non-identical current paths for both field directions as expected for a 
real disordered system. The voltage at the pair C-D at fields below ~ 0.3 T 
indicates that a longitudinal MR signal mixes into the overall signal showing 
an anomaly at the MIT at ~ 0.1 T [33]. 

We would like to note that an even Hall effect that occurs in inhomogeneous 
materials is a well known phenomenon. This effect has been observed in the 
mixed state of strongly inhomogeneous type-11 superconductors [38,39,40] as 
well as in inhomogeneous semiconductors [41,42]. 

In summary, transport measurements performed in different configurations on 
HOPG samples provide evidence for an integer-like QHE in agreement with 
theoretical expectations. Our experiments indicate that the absence of minima 
in the longitudinal MR (that usually accompany the plateaus in the Hall 
resistivity) as well as the non-saturating LMR observed in all graphite samples 
is due to the contribution of a Hall resistance generated by perpendicular 
current paths in the disordered media. This appears to be the origin for the 
new phenomenon reported here, namely the even QHE. 
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